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Abstract. 

We present a stellar populations analysis of the first release of the CFHTLS (Canada-France-Hawai Telescope 
Legacy Survey) data. A detailed comparison between the Besangon model of the Galaxy and the first data release 
of the CFHTLS-Deep survey is performed by implementing the MEGACAM photometric system in this model 
using stellar atmosphere model libraries. The reliability of the theoretical libraries to reproduce the observed 
colours in the MEGACAM system is investigated. The locations of various stellar species like subdwarfs, white 
dwarfs, late-type and brown dwarfs, binary systems are identified. The contamination of the stellar sample by 
quasars and compact galaxies is quantified using spectroscopic data from the VIMOS-VLT Deep Survey (VVDS) 
as a function of i' magnitude and r' — i' colour. A comparison between simulated counts using the standard IMF 
at low masses show that the number of very low mass dwarfs may have been underestimated in previous studies. 
These observations favour a power law IMF following d{n)/dm oc m~" with a = 2.5 for m < 0.25 M© or a = 3.0 
for m < 0.2 M0 for single stars. The resulting LF is in agreement with the local LF as measured from the 5 or 
25 pc samples. It is in strong disagreement with the Zheng et al (2001) LF measured from deep HST data. We 
show that this discrepancy can be understood as an indication of a different IMF at low masses at early epochs 
of the Galaxy compared to the local thin disc IMF. 

Key words. Galaxy: stellar content - stars; luminosity function, mass function -stars: white dwarfs - binaries: 
general -stars: stellar atmospheres -stars: low mass, brown dwarfs 



1. Introduction 

The CFHTLS (Canada-France Hawaii Telescope Legacy 
Survey) is a five year large observing program at the CFH 
Telescope, using the wide field prime focus MegaPrime 
equipped with MEGACAM, a 36 CCD mosaic camera. 
Together with its small pixel scale of 0.185 arcsec and the 
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large number of nights dedicated to the survey (around 
500 nights over five years), the CFHTLS goes deeper and 
has a better image quality than the Sloan Digitized Sky 
Survey but on a much smaller area of the sky. Hence it 
probes a different volume of the Universe. The scientific 
goals from CFHTLS cover a wide range of scales: from 
the solar system (systematic search of trans-Neptunian 
objects), stellar populations and galactic structure, up to 
the distant universe, constraining the geometry of the uni- 
verse (SNIa and cosmic shear), dark matter properties, 
quasars, clusters of galaxies, properties of galaxy cluster- 
ing and galaxy evolution at high redshift, and probing the 
relation between dark and luminous matter. 

We plan a series of papers dedicated to analysing stel- 
lar populations in the different surveys components of 
the CFHTLS. In this paper we analyse a subset of the 
first data release of the CFHTLS, investigating the ob- 
jects classified as stellar in the catalogues. We examine 
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Table 1. Location, median seeing and field of view (in sq. degree) of the 3 CFHTLS fields studied here. 



Field 


a(J2000) 


5(J2000) 


l{deg) 


b{deg) 


seeing in i' 


FOV 


Dl 


02h26m00s 


-04d30m00s 


172.0 


-58.0 


0.88 


0.80 


D2 


10h00m29s 


+02dl2m21s 


236.8 


42.1 


0.95 


0.69 


D3 


14hl9m28s 


+52d40m41s 


96.2 


59.6 


0.92 


0.77 



the photometric quahty of the three fields of the Deep 
Survey, Dl, D2 and D3, determine the stellar populations 
and the contamination by compact galaxies and quasars. 
Using the Besangon model of stellar population synthe- 
sis (Robin et al. l2U(J3p together with stellar atmosphere 
models BaselS.l (Lejeune et aL HWl Westera et al. l^TO^ 
and NextGen (AUard et al. ll997|l . we are able to produce 
synthetic star counts, colour-colour and colour-magnitude 
diagrams in the MEGACAM filter system. We determine 
the location of white dwarfs, brown dwarfs and binaries in 
the different colour-colour diagrams. We emphasize in this 
paper the study of the luminosity function at low masses 
of the disc population. A comparison between star counts 
of the M dwarf populations in the three CFHTLS Deep 
fields and model predictions provide new constraints on 
the slope of the IMF of low-mass stars, a parameter which 
is still under debate. 



2. The CFHTLS fields 

The CFHTLS survey consists of three different surveys: 

— The Deep, a survey comprising four fields (named Dl, 
D2, D3, D4) of 1 square degree each, in five filters 
u* ,g' ,r\i' and 2' reaching r' up to 28™, observed every 
3-4 days, 

— the Wide, a survey of 3 patches (Wl, W2, W3) each 
around 7 deg^ in five bands, with limiting magnitudes 
r' up to 25™, observed in two epochs separated by 3 
years, 

— the Very Wide, a survey of 1300 deg^ along the eclip- 
tic on 5-6 epochs, dedicated to the detection of fast 
moving solar system objects and stellar populations. 

We refer the reader to the corresponding web pages 
|http://www.c£ht. hawaii.edu/Science/CFHLS/) for a de- 
tailed description of the survey components. The observa- 
tions are carried out at the CFHT using the MEGACAM 
camera, which consists of a mosaic of 36 CCDs, each of 
them 2048 x 4612 pixels large. The pixel scale is 0.185" 
which gives a total field of view of 0.96 x 0.94 deg^. The 
seeing of the observations used to construct the stacks is 
better than 1.1", which guarantees high quality data for 
the three different surveys. 

We present in this paper an analysis of three fields 
of the "Deep" Survey, Dl, D2 and D3 of the CFHTLS 
release TOOOl. These data are stacks of many images. Field 
coordinates and the median seeing in i' in each field are 
given in Table 



3. Data reduction 

The stacks and catalogues used in this paper were released 
as part of the TERAPIX TOOOl pubhc release. A brief 
outline of how these stacks were prepared is as follows. 

CFHTLS observations are carried out with 
MEGACAM in queue survey mode. For release TOOOl, 
only observations from June, 1st 2003 to July, 22, 2004 
were used. Pre-reductions were carried out at the CFHT 
using the ELIXIR^ pre-reduction system at CFHT and 
then these pre-reduced images were shipped to TERAPIX 
via the Canadian Astronomy Data Centre, in Victoria, 
Canada. These pre-reduced images were then injected 
into the TERAPIX pipeline for inspection and quality 
control purposes. The TERAPIX tool QualityFITS was 
used to inspect and grade each image, and also to produce 
weight-maps derived from the CFHT-provided master 
fiats using the Weight Watcher tool. The global astromet- 
ric and photometric solutions were computed using the 
WIFIX package, an earlier generation of the TERAPIX 
astrometric software, as the production astrometric 
software package was still in testing phase at the time 
of the TOOOl release. For inclusion in the stacks, images 
must have a seeing better than 1.1" (1.3" in u*) and 
airmass less than 1.5. From this point on for the image 
reductions, we followed essentially the same procedure 
as outlined in McCracken et al. H2()()3|l . and we refer to 
the interested reader to this paper for more details. The 
two significant differences are firstly that we use weight 
maps computed from the image flat-fields themselves 
and secondly we use the USNO-B as the astrometric 
reference catalogue (which increases the robustness of the 
overall astrometric solution with respect to the solutions 
utilising the USNO-A). Full details of the properties of 
the final stacks, including depth in each filter and the 
accuracy of the final astrometric solution can be found 
on the TERAPIX web pages^,'^. The internal accuracy 
of the astrometric solution (band-to-band) is better than 
one pixel rms over the entire MEGACAM field, whereas 
the external astrometric solution is around ~ 0.25" rms. 

Photometric calibrations for each pre-reduced image is 
provided by the ELIXIR pipeline. ELIXIR also applies a 
CCD-to-CCD flux scaling derived from repeated observa- 
tions of dense stellar fields which are shifted many times 
around the MEGACAM field of view (providing magni- 
tude measurements of the same star on different CCDs). 
This procedure is necessary to correctly account for the 

^ http:/ /www. cfht.hawaii.edu/Instruments/Elixir/ 

http:/ /terapix.iap.fr/article.php?idarticle=383 
^ http:/ /terapix.iap.fr/article.php?idarticle=382 
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"scattered light" effect and ensures that the flux of any 
given object is independent of the position on the mo- 
saic. The residual ccd-to-ccd magnitude error following 
this procedure is around ~ 3%. In constructing the final 
stacks, we compare the magnitudes of objects in overlap- 
ping pointings and in each band the photometric expo- 
sures are indentified as those in which the objects have 
the highest flux: other images are scaled to these observa- 
tions. Based on an examination of galaxy counts and stel- 
lar colour-colour plots (see below), we estimate that our 
absolute photometric solution in each filter has a maxi- 
mum systematic error of ~ 0.05 magnitudes. Catalogues 
were extracted using SExtractor in dual-image mode, with 
detections carried out using a chi-squarcd image (Szalay 
et al. I2008|l constructed from the g'r'i' images. Kron-like 
total magnitudes were used throughout. Through this pa- 
per, our magnitudes are presented in the MEGACAM in- 
strumental AB system. 

4. Star-Galaxy separation 

We separated point-like sources from extended ones using 
SExtractor's (Bertin & Arnouts I1996| "flux-radius" pa- 
rameter measured on the z'-band image. This parameter 
measures the radius which encloses 50% of the object's 
flux: for point-like sources this is independent of magni- 
tude, and depends only on the image FWHM. The stars 
were selected by automatically locating the stellar branch 
in the flux-radius-magnitude diagram in a series of 10 arc- 
minute cells distributed over each MEGACAM i'-stack, 
which accounts for variation of FWHM over MEGACAM 
field of view. Figure ^ shows the compactness parameter 
against the magnitude for the three CFHTLS fields Dl, 
D2 and D3. At magnitudes fainter than i' = 21.0 the sepa- 
ration between stars and galaxies starts to be problematic. 
From Fig. n] it is clear that for the Dl field the contami- 
nation of galaxies is small for i' < 21.0 while for the D2 
and the D3 field the star/galaxy separation starts to fail 
at already i' ~ 20.5. The choice of the cut-off at i' — 21.0 
for the TOOOl release is certainly a conservative criterion 
which can not be applied for all three fields. However, the 
star/galaxy separation depends on the colour of the ob- 
jects. If one restricts to red objects with r' — i' > 1.4, stars 
can be better seperated from galaxies and thus stars can 
be extracted until i' < 22.5 (for the Dl field). We discuss 
and quantify below galaxy contamination as a function of 
magnitude and colour. 

In order to assess the number of galaxies contam- 
inating our stellar sample, we used spectroscopic data 
from the VIMO S-VLT Deep Survey (VVDS). The VVDS 
(Le Fevre et al. I2005|l intends to measure redshifts over 
< z < 5 across 16 deg^ in four separate fields. The sur- 
vey is conducted in two steps : an imaging survey and a 
spectroscopic survey. Deep imaging in the UBVRI bands 
and partly K' band is obtained with the CFHT, ESO- 
NTT and ESO-2.2m. The VLT VIMOS instruments al- 
low the measurement of redshifts of objects selected from 
the imaging survey. The so-called "Deep" survey has a 
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Fig. 1. i' magnitude vs half-light radius (in pixels) for the 
Dl, D2 and D3 field. The filled circles show the selected 
stellar branch. Note that objects with i' < 17 are satu- 
rated. 



limiting magnitude Iab = 24. Its location, a = 2 h 26 
m and S — — 4°30' overlaps with that of the Dl field. 
For a detailed description of the VVDs data we refer to 
Le Fevre et al. H2UU5f) . We extended the limiting mag- 
nitude for the star/galaxy separation of the Dl field to 
i' = 22.0 mag (thus Imag deeper than the official TOOOl 
release), using the half light radius, as described above. 
The Dl stacks have better seeing than the other fields, 
allowing star-galaxy separation to fainter magnitudes. We 
use this deeper catalogue in order to discuss the galaxy 
contamination (see below). However, for the rest of the 
paper we use the official TOOOl release, in which stars are 
separated from galaxies only until i' = 21.0. 

We cross-identified the Dl data with the VVDS data 
using a search radius of 2". Due to the masking and only 
partial overlap between the Dl field and the VVDS F02- 
field, only an area of 0.4 sq degree is in common. Out of 
9088 sources of VVDS data, 7110 sources have been cross- 
identified with the Dl field. 295 sources were identified as 
stellar and are shown in Fig. |2| 

Figure 121 shows the colour-colour diagrams of the Dl 
field cross-identified with the VVDS data. All objects 
shown are classified as stars from the morphological crite- 
rion. The dots show the spectroscopically identified stars 
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Table 2. Galaxy contamination of the Dl field cross- 
identified with the VVDS field as a function of i' mag- 
nitude and r' — i' colour. The common area is ~ 0.4 sq 
degree. 



r' -i' 


i' < 20.0 


i' < 21.0 


i' < 22.0 


total 


3.2% ±1.8 


6.7% ±1.9 


13.1% ±2.3 


< 0.5 


3.2% ±3.2 


13.0% ±4.3 


23.5% ±4.9 


> 0.5 


3.1% ±2.2 


2.7% ±1.5 


6.5% ±2.1 


> 1.4 


0% 


0% 


1.8% ±1.8 



to redder colour diminishes significantly the contamina- 
tion by galaxies. For r' — i' < 0.5 and i' < 22 one obtains 
for the Dl field a rate of 23.5% of contaminating galaxies 
while for r' ~ i' > 1.4 the percentage of galaxy contam- 
ination is negligible. Note that the galaxy contamination 
is larger for the D2 and D3 field. 

In addition, we visually inspected (for the Dl field) 
sources which are located outside the stellar locus. We 
have marked them in Fig.|21as open squares. Most of these 
objects are galaxies blended with stars where obviously 
the aperture photometry is unreliable. Note that only the 
g' — r' vs r' — i' diagram reveals these blends easily. 
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Fig. 2. Colour-colour diagrams of the Dl field cross- 
identified with the VVDS data. Objects shown are all mor- 
phologically classified as stars. Dots are spectroscopically 
identified stars. Asterisks are galaxies and open squares 
indicate sources where a blend of a galaxy with a star has 
been found by visual inspection of the images. The lim- 
iting magnitude used for the morphological star/galaxy 
separation is i' = 22.0 mag. 



and asterisks those objects spectroscopically identified as 
galaxies. The majority of the galaxies classified morpho- 
logically as stars populate the blue part of the colour- 
colour diagram with r' — i' < 0.4. They populate the stel- 
lar locus as well as regions outside of the stellar sequence. 
Table [21 gives the percentage of galaxies contaminating 
the stellar sample as a function of i' magnitude and r' — i' 
colour. The percentage of galaxy contamination depends 
very much on the r' — i' colour, which means that going 



5. Stellar atmospheres and the theoretical stellar 
locus 

Using the published response of the CFHT, Megaprime 
and MEGACAM we have computed synthetic colours of 
the stars as a function of temperature, gravity and metal- 
licity. To do this we have used two sets of stellar atmo- 
sphere models: Basel 3.1 and NextGen. The Basel3.1 li- 
brary is a semi-empirical library based on preceding gen- 
eration of models, Basel 2.2 (Lejeune et al. I1997fl . ex- 
tended to non-solar metallicities by Westera et al. H2002|l . 
The Kurucz theoretical spectra H1979|l have been modi- 
fied to fit broad band photometry using the algorithm by 
Cuisinier et al. (see Buser& Kurucz [l992|l . The corrected 
spectra are used for integrating the flirx in the desired 
bands. NextGen is the 1997 version of atmosphere models 

mm 



from AUard et al. H1997|l . These models use a direct opac- 
ity sampling including over 500 million lines of atomic and 
molecular species. They give a more realistic description 
of the M dwarf population. 

Figure |21 shows the CFHTLS colour-colour diagrams 
of the Dl field superimposed with the synthetic colours 
of dwarf stars using the Basel3.1 stellar library for solar 
metaUicity, [Fe/H]= -1.0 and [Fe/H]= -2.0 (right panel) 
and the NextGen library for [Fe/H]=0.0 and [Fe/H]= -1.0 
(left panel). For a better definition of the stellar locus, we 
use here only stellar objects which have a photometric 
error estimate smaller than 0.01 mag in each filter. We 
noted that the i' — z' colour has a slight offset of 0.05 mag 
in the Dl field compared to the model which comes from 
uncertainties on the photometric calibration. For the D2 
and D3 field we noted a shift of 0.07 mag and 0.02 mag in 
r' — i' respectively. These offsets have been applied. 

Figure 3 illustrates the sensitivity of the colours in the 
CFHTLS system to metallicity and the differences in the 
stellar libraries. 

For temperatures below 3500 K, which correspond to 
K/M stars, the BaselS.l library does not give realistic 
colours for cool dwarfs, but gives a better fit to the data 
than NexGen models do for hotter stars. 

In the temperature range 7000 to 4000 K, the most 
sensitive colour is g' — r', whereas for cooler stars this 



colour index saturates and 



becomes better. The 



colour seems to be redundant with r' —i' , but going 




Fig. 3. Colour-colour diagrams of the Dl field compared to the synthetic colour-colour diagram of NextGen (left panel) 
and BaselS.l (right panel) for different metallicities. The solid line denotes [Fe/H]=0.0, the dashed line [Fe/H]= —1.0 
and the dotted line [Fe/H]= —2.0. For a better determination of the stellar locus, we show here only stellar objects 
with a photometric error smaller than 0.01 mag in each filter. The approximate temperatures are also indicated. 



to very cool stars we expect it to be a very good indicator 
for selecting brown dwarfs (see sect.7.2). 

Both atmosphere models show a strong metallicity ef- 
fect for cool dwarfs in g' — r' and i' — z' . It appears that 



at 5' — r' > 1 this index is no longer sensitive to temper- 
ature but to metallicity. If the photometric calibration of 
the survey is accurate and the model atmospheres reliable, 
we will be able to constrain the metallicity distribution of 
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model 




Fig. 4. Colour-colour diagram predicted by the Besangon 
model for the Dl field for thin disc (dots), thick disc (open 
squares) and spheroid (asterisks) stars. 



these cool stars, at least statistically. This will permit us 
to determine the metallicity range, and probably the thin 
disc to thick disc density ratio, as we expect a difference 
of metallicity of about 0.5 dex between these two popula- 
tions, corresponding to about 0.15 magnitude in g' — r' at 
a temperature of 4000 K. 

6. The Besan^on Galaxy model 

The Besangon Galaxy model is a simulation tool aimed at 
testing galaxy evolution scenarii by comparing stellar dis- 
tributions predicted by these scenarii with observations, 
such as photometric star counts and kinematics. A com- 
plete description of the model ingredients can be found 
in Robin et al. H2003|l . We summarise here the model's 
principal features. 

The model assumes that stars are created from gas fol- 
lowing a star formation history and an initial mass func- 
tion; stellar evolution follows evolutionary tracks. To re- 
produce the overall galaxy formation and evolution we 
distinguish four populations of different ages and star for- 
mation history, which we now describe. 

The oldest population is the spheroid. For this popula- 
tion we assume a single-burst star formation history ocur- 



ring early in the lifetime of the Galaxy, around 14 Gyr 
ago, from gas still in a spheroidal distribution. The initial 
mass function (IMF) and the density distribution of this 
population is characterised by a power law index, an axis 
ratio and a local normalisation, which are constrained by 
remote star counts at high and medium Galactic latitudes 
(Robin et al. l2UUU|l . The kinematics are also deduced from 
in situ velocity determinations. The mean metallicity of 
the spheroid is assumed to be -1.7dex with an intrinsic 
dispersion of 0.5 dex. 

Secondly, a population, called the thick disc, is formed 
of stars born about 11-12 Gyr ago in a short period of time 
as implied by recent metallicity determinations for this 
population. We also assume a single burst for simplicity. 
For the thick disk, star formation occurred from the gas 
already settled in the disk. The kinematics, deduced from 
observational constraints (Ojha et al. 119961 Oiha ITM^ . 
implies that it has undergone a merging event shortly after 
the disc formation (Robin et al. I1996|l , increasing the disk 
thickness and giving a higher velocity dispersion and scale 
height. The IMF, density distribution and local normali- 
sation were constrained from star counts (Reyle & Robin 
12001(1 . The mean metallicity of the thick disc is assumed 
to be -0.7 dex with an intrinsic dispersion of 0.3 dex. 

Thirdly, a bulge population is present in the center 
of the Galaxy and extends to about 2 kpc. Its age is of 
the order of 10 Gyr. This value is however poorly con- 
strained. This population has a triaxial distribution, as 
a bar as determined by Picaud et al. f2004') from near- 
infrared star counts. Velocity dispersions are large, simi- 
lar to the spheroid. The mean metallicity is assumed solar 
with a dispersion of 0.5 dex. 

A standard evolution model is used to compute the 
disc population, based on a typical set of parameters: an 
IMF, a star formation rate (SFR), a set of evolutionary 
tracks (see Haywood et al., 1997a and references therein). 
The disc population is assumed to evolve during 10 Gyr. A 
set of IMF slopes and SFR's are tentatively assumed and 
tested against star counts. The tuning of disc parameters 
against relevant observational data has been described in 
Haywood et al. H1997al[T997Ejl . 

A revised IMF has been used as a starting point in the 
present analysis, adjusted to agree with the most recent 
Hipparcos results: the age- velocity dispersion relation is 
from Gomez et al. I)1997|l . the local luminosity function 

IM3 



is from Jahreifi & Wielen H1997|l and an IMF is adjusted 
to it, giving an IMF slope a = 1.5 in the low mass range 
[0.5-0.08 Mq], in good agreement with Kroupa (,,2001,) . The 
scale height has been computed self-consistently using the 
potential via the Boltzmann equation. The local dynami- 

mm. 



cal mass was taken from Creze et al. 

The evolutionary model fixes the distribution of stars 
within the parameter space of effective temperature, grav- 
ity, absolute magnitude, mass and age. These parameters 
are converted into colours in various systems through stel- 
lar atmosphere models. 

In the case of the MEGACAM photometric system, 
we have used the optics, CCD and filter definition of the 
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passbands, and applied them to the spectral libraries. As 
seen in Fig. 1, the BaselS.l library is more suitable for hot 
stars, giving better predicted colours, especially u* — g' , 
while the NextGen library is more realistic for cool stars. 
Hence we have adopted a combination of both: BaselS.l 
at Teff > 4000 K and NextGen for cooler stars. As the 
cooler stars are mostly located in the disc, we used only 
[Fe/H]=0.0 and [Fe/H]= -1.0 for Tos < 4000 K. The 
model simulations also include a model of extinction and 
account for observational errors. The Besangon model pre- 
dictions in the MEGACAM photometric system can be 
found at jhttp: / /www.obs-besancon.fr /modelCj 

7. Stellar populations 

Simulations from the Besangon model in the MEGACAM 
photometric system have been performed and compared 
with the CFHTLS in three fields: Dl, D2 and D3. 

We used the same magnitude limit in i' of 21.0 im- 
posed by the star/galaxy separation and took the pho- 
tometric errors in each band as a function of magnitude 
into account. In the discussion below we indicate also the 
different components of the Galactic model such as thin 
disc, thick disc and the spheroid. Figure^ shows the ex- 
pected colour-colour diagram predicted by the Besangon 
model for the Dl field. Indicated are the three different 
components, thin disc, thick disc and halo. While the 
blue part of the colour-colour diagram is populated by 
mainly spheroid stars, the thick disc stars are concentrated 
around 0.6 < r' — i' < 0.9 and the thin disc stars cover 
the red part of the diagram. 

Figure El shows the colour-colour diagrams of the 
CFHTLS for the Dl, D2 and D3 fields respectively as weU 
as the model predictions. Note the excellent overall agree- 
ment between observed and predicted colours for the three 
Galactic components. Figure [Sjalso shows that the stellar 
populations of the three deep fields are similar, although 
one notices several differences: 

— there is a steady increase in the width of the stellar lo- 
cus from the Dl field - which has a well defined narrow 
locus - towards the D3 field - which is wider - and the 
D2 field which is wider still, especially at r' — i' > 1.2, 
and more pronouced in i' — z' . 

— The blue part of the colour-colour diagrams {g' — r' < 
0.5, r' — i' < 0.4, i' — z' < 0.3) is populated by an 
ovcrdcnsity of objects which has a larger spread than 
the predicted stellar locus. These objects are compact 
faint galaxies or quasars erroneously classified as stars 
by morphological criteria (see Sect. 4). This contami- 
nation is larger in D2 and D3 than in Dl. 

— In the g' — r' vs r' — i' diagram, there are sources with 
g' — r' > 1.3 and r' — i' < 1.5 abundant in the D2 
and D3 fields, but absent in the Dl field. Pello (priv. 
communication) provided us with theoretical colours 
of galaxies in the MEGACAM photometric system for 
different redshifts and for different morphological types 
which show that these objects are expected to be el- 



liptical galaxies with z < 0.5 while a few high redshift 
QSOs and spiral galaxies could also be present. 

— In the Dl, D2 and D3 fields there are objects with 
r' — i' > 2 and i' — z' > 1.2 which are off the stellar 
sequence which may be brown dwarf candidates (see 
Sect. 7.2), or high redshift quasars. 

— The dispersion in the synthetic colour-colour diagram 
is larger than observed, especially in the blue part. 
This is due to uncertainties of the model atmospheres 
in this metallicity range, as seen in Fig. 3. 

7.1. White dwarfs 

White dwarfs (WD) are the last stage of stellar evolu- 
tion and their space density depends on the Galactic star 
formation history and initial mass function. While the lu- 
minosity function (thereafter LF) of white dwarfs in the 
thin disc is known from systematic searches in the solar 
neighbourhood (e. g. Liebert et al. 119881 Ruiz & Bergeron 
2001, Holberg et al. 12002(1 . only very few thick disc white 
dwarfs have been identified while the presence of white 
dwarfs in the Galactic halo is still uncertain. Knowledge 
of the luminosity function of thick disc and halo white 
dwarfs is expected to constrain the age of these popula- 
tions, the physics of the coolest white dwarfs, as well as 
the initial mass function at early epochs in the Galaxy 
(ChabrierEOnSll. 

Creze et al. H2004|l used two-epoch observations of the 
1 deg^ VVDS-F02 deep field to search for white dwarfs in 
the VVDS survey by proper motions; they reported a null 
detection. 

Bergeron (priv. communication) provided us with the- 
oretical colours of white dwarfs in the MEGACAM pho- 
tometric system. Figure |S1 shows their location in vari- 
ous colour-colour diagrams compared with data in the Dl 
field. The locus of the white dwarfs is distinct from the 
locus of the subdwarfs in the u* — g' vs g' — r' plane only. 
Furthermore, we note the sensitivity of the colours to logg. 
However this part of colour-colour space is highly con- 
taminated by quasars. In the other colour combinations, 
white dwarfs are indistinguishable from subdwarfs, except 
the very cool ones (temperature less than about 3200 K) 
where DA white dwarfs with hydrogen atmospheres start 
to deviate from blackbody. 

In the present data, there are a few objects with colours 
consistent with those of white dwarfs. However these ob- 
jects could also be horizontal branch spheroid stars or 
quasars. White dwarf candidates will only be reliably iden- 
tified when proper motions become available. Proper mo- 
tions will easily distinguish these objects from horizontal 
branch stars, as they are brighter, hence much more dis- 
tant at a given apparent magnitude (extragalactic objects 
such as quasars, of course, will have no measurable proper 
motions). The CFHTLS will allow us to eventually cover 
about 150 square degrees to the same magnitude limit 
with proper motions, allowing a definitive answer if bary- 
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Fig. 5. Colour-colour diagrams of the three Deep field Dl, D2 and D3 compared with synthetic colours predicted by 
the Besangon model (lower panel). Note the excellent agreement between the data and the model. 



onic matter is present in galactic hales in the form of white 
dwarfs. 

The number of expected white dwarfs per square de- 
gree, as predicted by the Besangon Galaxy model is about 
25 for the thin disc, around two for the thick disc and 



about 0.1 for the halo, to magnitude i' = 22.5. The num- 
ber of ancient halo WDs is computed assuming that the 
dark halo is partly made of ancient white dwarfs, at the 
level of 2% of its mass density. Even at these bright mag- 
nitudes their identification from photometry only will be 
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Fig. 6. White dwarf sequences in colour-colour diagrams 
of the Dl field. Solid line shows the DA sequence, 
dashed line the DB sequence from theoretical model from 
Bergeron et al. Ij^UUlfl for surface gravities of logg=7.0, 
8.0 and 9.0 



difficult due to the large number of galaxies and quasars; 
proper motions are necessary. Assuming that the survey 
would reach an astrometric accuracy of about 0.1 pixel at 
i' = 22.5, with a time baseline of three years, objects hav- 
ing a motion of 25 mas/yr would be detectable at the S-cr 
level. About 2/3 of the disc WD pass this proper motion 
selection criterion, 80% of the thick disc WDs, and 100% 
of the halo WDs. 



7.2. M, L, and T dwarfs 

Large sky surveys are very efficient in identifying ex- 
tremely late-type stellar and sub-stellar objects because 
they are detectable only at short distances from the sun, 
and hence in a small volume. Even if they are numerous, 
the number of brown dwarfs detectable to date is still 
small. Near-infrared surveys such as DENIS and 2MASS 
resulted in the detection of a significant number of close 
low-mass stars and brown dwarfs (see e. g. Kendall et 
al. 'lUUM Reid et al. VMM Burgasser et al. OUMjl . The 
SDSS survey was used to identify the first field T dwarfs 
(Strauss et al. 119991 Tsvetanov et al. I2000|l and the num- 
bers of known L dwarfs has been greatly increascd(Fan et 
al. nUOO. Schneider et al. lMI^ . Gelino et al. ( 2004j intro- 
duced an homogeneous database of M, L, and T dwarfs, 
that contains more than one thousand of these objects. 
Hawley et al. ( 2002) compiled a large sample of M, L and 
T dwarfs from SDSS spectra together with SDSS photom- 
etry and additional near-IR photometry (2MASS). They 
find that the i — z and i — J colour are the most useful 
for estimating spectral types based solely on photometric 
information for M and L dwarfs. 

Figure[71shows the i' — z' vs r' — i' diagram of the three 
CFHTLS fields together superimposed with the NextGen 
model at solar metallicity. The spectral types are from the 
temperature-spectral type relation from Golimowski et al. 
( 2004 ) . There are two objects with r'-i' > 2 and i' - z' > 
1.3 in the D3 field and one candidate each in the Dl and 
D2 fields, which arc good candidates for being either early 
L dwarfs or high redshift quasars (see Fig. d). They are 
given in Tab. 3. High redshift quasars with i' — z' > 1.5 
will be distinguishable from brown dwarfs either by near- 
infrared photometry or by proper motion measurements. 



7.3. Binary systems 

The formation and evolution of low-mass stars in a binary 
system is a common phenomenon which leads to the in- 
teresting class of cataclysmic variables. In deep surveys 
one expects to detect a few cases of WD-M dwarf pairs. 
Raymond et al. P003) identified ~ 100 white dwarf-M 
dwarf pairs in the SDSS survey with g < 20. Using addi- 
tional spectroscopy, they achieve an efficiency of ~ 60% 
in finding white dwarf-M dwarf pairs because of the con- 
tamination by galaxies in the interesting colour regions. 

We simulated a sample of unresolved M dwarfs + white 
dwarfs systems by merging their fluxes. Typical colours of 
these simulated systems are given in Fig.|Slas star symbols. 
The location of these systems is clearly outside the single 
star locus in the g' — r' vs r' — i' diagram. However they 
lie in a region where we expect contamination by compact 
galaxies and quasars. Their identification will be easy us- 
ing proper motions, all these objects being intrinsically 
faint, and are hence detected only in the solar neighbour- 
hood. 
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Table 3. Positions and MEGACAM photometry of late type dwarfs or high redshift quasars. 
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Fig. 7. r' — i' vs i' ~ z' diagrams of the Dl, 
fields. Superimposed are the average colors c 
dwarfs taken from NextGen models. Spectral t 
dicated, following the temperature-spectral tj 
from Golimowski et al. |20Q4| ). Candidate L 
shown as asterisks (see also Tab. 3). 

8. Stellar densities and the IMF at low r 

As mentioned above, the stellar samples are co 

by non-stellar sources. As a large fraction of t' 

side the stellar locus and stellar binaries and w 

are also expected outside the main sequence ( 

only proper motions can be used to clean our 

remove galaxies and quasars. In the following all objects 

classified as stellar are kept. 

In Fig. |51 histograms in g' — r' , r' — i' and i' — z' for 
the Dl field are shown, with model predictions for each 
population, thin disc, thick disc, and spheroid. Model pre- 
dictions are acceptable for all three populations, which are 
better separated in the r' — i' colour. 

8.1. Stellar densities 

Figure 1101 shows the g' — r' vs r' — i' diagram for each 
CFHTLS field, compared with model predictions, where 
the greyscale indicates the number of stars. The colour- 
colour diagrams in the three fields are similar, and model 
predictions are in good general agreement with the data. 
We notice, however, a few significant differences: 
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Fig. 8. Binary systems M dwarf -|- White dwarf in the D3 
field. Dots: whole Dl data set. Stars: simulated systems. 
Colours are estimated by adding up the flux in a realistic 
sample of unresolved systems. 



— In all three fields the predicted source density is a 
slightly low compared with the observations. The dif- 
ference between the data and the predicted colour 
counts is of order ^ 30% for r' — i' < 0.3. We estimate 
that part of the excess in the data is overestimated, 
due to the contamination by galaxies, which in the 
blue is of order ^ 13% for the Dl field but larger for 
the D2 and D3 field (see also Section 4). On the other 
hand, the predicted counts of spheroid stars depend on 
three density parameters: the axis ratio, the power law 
exponent and the local normalization. These parame- 
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Fig. 10. g' — r' vs r' — i' diagram of the Dl, D2 and D3 



field compared to the synthetic colours predicted by the 
Besangon model, all populations included. The greyscale 
indicates the number of stars per bin width of 0.05 mag in 
Fig. 9. Histograms of the colour distributions in the Dl each colour, 
field. The thick hne denotes the observations, the thin line 
the Besangon model, the dotted line the contribution of 
the spheroid, the short-dashed line the thick disc and the 
dot-long dashed line the thin disc contribution. 



ity function of these stars has been underestimated in 
the model. We investigate this point further below. 



ters have been already constrained by other data sets 
(Robin et al. l2000|l but could be better adjusted using 
new large and deep surveys such as the CFHTLS or 
SDSS. In order to perform better adjustments of these 
parameters, we need more fields with wider areas of the 
sky than those presented here. The effects on the star 
counts of the derived IMF of this population as well 
as their density distribution can only be disentangled 
from large surveys in different areas of the sky. This 
will be considered in a future paper using the T0002 
and T0003 releases. 

— In the range 0.3 < r' — z' < 0.7 the thick disc popu- 
lation dominates while for r' — i' > 0.8 the thin disc 
dominates the counts. The whole colour range is well 
modeled in the three fields; the thick disc model fitted 
in Reyle & Robin (2001) fits well here too, in addition 
to the thin disc for 0.8 < r' - i' < 1.4. 

— In the very red part of the diagram, at r' — i' > 1.6, the 
number of observed stars is significantly larger than 
predicted by the model. We do not expect large num- 
bers of very high redshift galaxies and quasars at these 
colours. These point source objects are most probably 
disc M dwarfs. This implies that the assumed luminos- 



8.2. The IMF at low masses 

The behaviour of the stellar luminosity function (LF) and 
mass function (MF) for low-mass stars (< 1 Mq) is still 
under debate. Jahrcifi et al. H1997|l derived a local stellar 
LF (within 20 pc) from the Catalogue of Nearby Stars 
revised with Hipparcos data. Reid et al. (lU^) derived 
the nearby star LF using the Palomar/Michigan State 
University sample (PMSU) combined with the Hipparcos 
25 pc sample. Both LFs are shown in Figure lTTl Error bars 
are large due to the small survey volumes. For masses 
smaller than 0.5 M0, the determination of the MF is 
hampered by the incompleteness of the different samples 
(Henry et al.dSS?, Chabrier & Baraffe 2000:) and by the 
unknown proportion of M dwarfs in binaries. Chabrier 
1)200311 estimated that the mass function below 1 Mq is 
consistent with a fraction of ^ 50% of M dwarf binaries 
where 30 % should have an M dwarf companion and 20% 
a brown dwarf secondary. 

Also on Figure 1111 are superimposed the luminosity 
functions used for simulations with our standard Galaxy 
model, as well as a few other luminosity functions obtained 
by varying the IMF slope at low mass. The luminosity 
function is made from segments of a power law IMF, as 
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luminosity function strongly depends on the assumed IMF 
slope and on the mass at which the slope changes: 



Fig. 11. Luminosity function in the V band in number of 
stars per pc^ per magnitude for different break points in 
mass. The top panel indicates nic = 0.15 Mq, the middle 
panel rric — 0.20 Mq and the lower panel rric = 0.25 M0. 
Open squares are from Jahreiss et al.'s H1997|l determi- 
nation from the revised Catalogue of Nearby Stars, open 
circles are from Reid et al. H2002|l determination using 
the PMSU survey combined with Hipparcos data, filled 
squares are from Zheng et al.'s ( 20011 determination from 
the HST. The lines show model luminosity functions as- 
suming different slopes for the IMF: a = 1.5 is the stan- 
dard Galaxy model (dashed line), a — 2 (thin line), a = 3 
(thick line), and a = 2 (dotted line) (see text). 



given in equation Q and a mass luminosity relation from 
Delfosse et al. (|2000)) in the magnitude range 12< A/y < 
17 and from theoretical models of Baraffe et al. (1998) 
at lower masses. The cutoff between absolute magnitude 
16 to 18 is mostly due to the mass-luminosity relation 
and only weakly dependent on the IMF slope at the very 
low mass end. However in the range 13< My < 16 the 



dn 

- — oc m 

am 

with: 

a = 1.5forTO<0.5 Mq (standard Galaxy model) 
a — 1.5 for m < 0.5 Mq and a — 2 for m < rric 
a — 1.5 for m < 0.5 Mq and a = 3 for m < rric 
a = 1.5 for m < 0.5 Mq and a = 4 for m < mc 



(1) 



and rrir = 0.15, 0.20, 0.25 Mq 



0- 



In the following the various tested IMFs are denoted 



LFa 



where alpha is the IMF slope , and rric is the 



mass where the slope changes. 



We note that simulated stars considered here are sin- 
gle stars. Distant binary systems may be not correctly 
identified as stars, although from a detailed analysis of 
the binary effect (a complete analysis of the binary effect 
is postponed to the next paper of this series), we have 
estimated that at r' — i' > 1.6 (that is My > 14) the 
correction for stars missed for this reason is negligible. At 
r' — i' < 1.5 the correction would be less than 25% and at 
1.5 < r' — i' < 1.6 it is less than 13%. Hence with respect 
to the very faint end of the LF studied here, the binary 
effect is expected to be negligible because these stars are 
too close to have been missed even in binary systems. We 
limit our further comparison to r' — i' > 1.5. 

Figure [TUI shows clearly a deficit of late type thin disc 
dwarfs in model predictions at r' — i' > 1.6, for the stan- 
dard Galaxy model. We have attempted to fit the LF to 
the available CFHTLS data in the 3 fields using the IMF 
formula from eq. (1). Figure IT^ shows the difference of 
the predicted star counts compared to observations for 
the three tested IMFs: a = 2 , a = 3 and a = 4 with 
mc = 0.20 Mq. 

In Table 4, 5 and 6 we give the number of stars in r' — i' 
intervals from the Dl, D2 and D3 field respectively, and, 
for comparison, the numbers simulated from tested LFs, 
varying TOc from 0.15 to 0.25 and a at m < mc from 1.5 
to 4. The standard deviation of the models relative to the 
data are estimated from Poisson statistics: assuming that 
the main source of error in the data is Poisson noise, we 
estimate the deviation of the model by computing : 



^model 



Ndata — 

V Ndata 



which gives the relative difference in the counts in units 
of the Poissonian scatter. Acceptable models have Umodei 
in the range [—3, 3]. 

From these Tables we can conclude that : 
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Table 4. Number of stars in the Dl field for i' < 21 and in different intervals of r' — i'. Column 1 indicates the middle 
of the interval of width 0.2, "tot" means the total of the 3 interval considered, that is 1.5 < r' — i' < 2.1. Columns 2 to 
17 give the model counts for each of the tested LF while column 18 contains the observed counts. Models are described 
by two parameters: on the first line of the column heading the value of the mass of changing slope is indicated, the 
second line gives the value of a (see eq^ . The column values include, in the first line, the number of stars in the colour 
interval, the second line gives the difference between the data and simulated counts in number of sigmas, assuming 
that the noise is dominated by the Poisson statistics (see text). Model with values between -3 and +3 are considered 
as acceptable in the colour bin. 
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Table 5. Same as Table lO but for field D2. 
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Table 6. Same as table El but for field D3. 



— The standard LF (with a = 1.5, i.e. no change of slope) 
is rejected at more than 6 sigmas in the three fields. 

— The LF with a rric = 0.15 Mq is deficient in stars 
compared with the data in all three fields, regardless 
of the slope. 

— Consistency with the data can be obtained with nic — 
0.25 or rric — 0.20 Mq. Depending on which to^ we 
take, the value of a is slightly changed. Having a higher 
rric allows a smaller a. Constraints from the local LF 
(see Figure 11) also indicate that having rric > 0.25 
Mq would not allow large values for a. 

— The counts in all 3 fields are best fitted by an IMF 
slope of a = 3 (with rric = 0.20 Mq) or a = 2.5 (with 



rric = 0.25 Mq). As seen from Figure 11 this LF is in 
acceptable agreement with the local LF as determined 
from the 25 pc sample. Comparing this IMF (either 
with mc = 0.20 or 0.25 Mq) with Chabrier ^^UU^ 
log-normal LF, yields good agreement but the revised 
IMF from Chabrier l|2004|l gives a too few low mass 
stars for 0.10-0.15 Mq compared with our determina- 
tion and with the local IMF. We emphasize that we 
have confirmed, from our intermediate distance sam- 
ple that the disc IMF does not decrease for masses 
above the 0.1 Mq limit. 
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Fig. 12. Difference between model and data as a function 
of the t' — i colour for the Dl, D2 and D3 field. The solid 
line indicates the standard model [a = 1.5) with mc = 
0.20 Mq, the dotted line assuming a — 2, the short dashed 
and the long dashed line assuming a = 3 and a = 4 
respectively, as explained in the text. The error bars are 
calculated assuming that the noise is due to the Poisson 
statistics. 



It should be noted that the number of stars decreases 
rapidly with the increasing r' — i' colour. A slight change 
in the zero point of the photometric calibration in r' or 
i' may change slightly our conclusions. We estimate from 
Tables 4, 5 and 6 that a systematic shift of 0.05 magnitude 
in r' — i' would produce a change of about 0.5 on the slope 
a. However such a shift is improbable as it would be seen 
also at the blue end of the sequence (r' — i' = 0) which is 
not the case, as seen in Figure 9. 



8.3. Comparison with HST results 

mm 



Zheng et al. H2001|l determined the luminosity function 
from a sample of about 1400 M dwarfs in 148 fields using 
the WFC2 and 162 fields from PCI with the HST. Their 
sample is characterized by a mean height above the plane 
of 1.5 kpc, with very few stars at vertical height <1 kpc. 
They derive their LF and IMF taking into account a prob- 
able metallicity gradient, by adopting a metallicity of -0.5 
dex at 1.5 kpc, and a colour-absolute magnitude relation 



varying with metallicity. Hence their sample is dominated 
by what we usually call the thick disc population. They 
deduce an IMF slope of a = —0.10 or a = —0.47 with or 
without the metallicity gradient taken into account. 

The sample considered in this paper is significantly 
different from the HST sample, as it is dominated by stars 
at distances above the plane of 150 to 450 pc with a mean 
distance of 350 pc for stars at r' — i' = 1.6 and 210 pc 
for stars having r' — i' = 2.0. This has two consequences: 
1) the sample is less biased by unresolved binaries and 2) 
it is dominated by the normal thin disc population and 
more comparable with the local sample which is used to 
determine the LF in the solar neighbourhood (Reid et al. 
1^)04.1 . 

Reyle & Robin f2001') have performed the first deter- 
mination of the thick disc IMF from a multi-directional 
analysis of star counts. They obtained an IMF dN/dm cx 
jj-^ ^j^g mass range 0.2 < m < 0.8Mq, which is in 
agreement with the IMF deduced by Zheng et al. H2001|l 
from the HST sample, reinforcing the ideas that: firstly 
Zheng et al. H2001I) have measured the thick disc IMF, 
rather than the thin disc one; secondly the thin disc and 
thick disc have different IMF slopes at low masses. The 
IMF found by Reyle & Robin l|MIT|l in the thick disc 
is well in agreement with the one determined in globular 
clusters (Paresece & De Marchi l2000f) and significantly dif- 
ferent from the one found in the local thin disc (Kroupa 
.20011 . 

The origin of the thick disc has long been a matter 
of debate. Nowadays favoured scenarii explain the thick 
disc by one or more accretions of galaxy satellites at early 
epochs of the Galaxy's formation, or by star formation 
from gas accreted during a chaotic period of hierarchical 
clustering (Brook et al. I2004|l . The thick disc is old and 
metal poor relative to the sun and it is also enhanced in 
alpha elements. Abundance determinations (Gratton et 
al. I2000|l also show that there has been a discontinuity 
in the star formation between the thick disc and the thin 
disc of at least 1 Gyr. The conditions of star formation 
at the epoch of thick disc formation were clearly differ- 
ent from the present conditions in the thin disc. Larson 
(2005) has analysed the physical conditions required for 
the thermal coupling of gas and dust in cloud fragmenta- 
tion. He studied the roles of the metallicity, background 
radiation and dust environment on the Jeans mass, hence 
on the typical mass of the stars formed. The combined 
effects of the metallicity and the possible lack of dust at 
cosniological epochs could increase the peak mass of the 
IMF relative to the present one. Moreover as the cosmic 
background temperature was higher in the past, a higher 
minimum cloud temperature exists, which also implies a 
higher Jeans mass. These conditions may well explain the 
fact that the IMF found in the thick disc has a typical 
mass higher than the thin disc, and is deficient in the very 
low mass stars which are found in the present disk MF. 
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9. Conclusions 

We have presented an analysis of the stellar populations in 
the CFHTLS using catalogues and images from the first 
public data release. Our population synthesis approach 
allowed us to test stellar libraries and to identify differ- 
ent stellar types and Galaxy components using colour- 
colour diagrams. We discuss the locations of various stellar 
species such as white dwarfs, late-type and brown dwarfs 
and binary systems in the MEG AC AM filter/detector 
combination. The contamination of the stellar sample by 
quasars and compact galaxies is quantified using spectro- 
scopic data from the VIMOS-VLT Deep Survey (VVDS). 
The percentage of the galaxy contamination depends very 
much on the r' — i' colour and can reach a maximum of 
23.5% for i' < 22.0 and r' - i' < 0.5. 

Our main conclusions concern the luminosity and mass 
function (MF) at low mass for the disc population. This 
data set favours an MF slope of a = 2.5± 1.0 for m < 0.25 
Mq or a = 3.0 ± 1.0 for m <0.2 Mq which although 
steep compared with previous investigations from other 
deep imaging surveys (such as Zheng et al. who used HST 
images) is still in agreement with local determinations of 
the IMF. This discrepancy can be explained by differences 
in the mean age and physical conditions of star formation 
of the samples, one being at about 1 kpc or more where the 
thick disc population is expected to dominate, and ours 
being at 150-450 pc and dominated by the thin disc. This 
discrepancy between the thin disc and thick disc IMFs 
could be explained if for physical reasons (for example, 
lack of dust, higher temperature backgound radiation or 
metallicity) very low mass star formation has been less 
efficient at the epoch of the thick disc formation. 

The new IMF as determined here cannot be extrap- 
olated to masses below 0.1 Mq. It is probable from the 
numbers of known brown dwarfs in clusters that the IMF 
starts to decrease near the H burning limit (Kroupa'200r) . 

In future papers we plan a more detailed analysis of 
these stellar populations, in particular the IMF at low 
masses of the disk, thick disk and spheroid and the old 
population density distribution up to several tens of kilo- 
parsecs. 

This might be performed using more accurate star- 
galaxy separation and by accounting for binary frequency 
in the modelling. The combination of multiband wide sur- 
vey coverage together with proper motions will enable us 
to count thick disc and halo white dwarfs, and to con- 
strain on the fraction of baryonic dark matter present in 
the form of stellar remnants. 
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